Topological Quantum Phase Transition in 5d Transition Metal Oxide Na 2 IrOa 
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We predict a quantum phase transition from normal to topological insulators in the 5d transition 
metal oxide Na2lr03, where the transition can be driven by the change of the long-range hopping 
and trigonal crystal field terms. From the first-principles-derived tight-binding Hamiltonian we de- 
termine the phase boundary through the parity analysis. In addition, our first-principles calculations 
for Na2lr03 model structures show that the interlayer distance can be an important parameter for 
the existence of a three-dimensional strong topological insulator phase. Na2lr03 is suggested to 
be a candidate material which can have both a nontrivial topology of bands and strong electron 
correlations. 

PACS numbers: 71.70.Ej, 73.20.-r, 73.43.Nq 



Topological insulators are newly discovered materials 
with a bulk band gap and topologically protected metallic 
surface states [U [2] ■ The theoretical predictions on Bi- 
based topological insulators, such as Bi x Sbi_ x , Bi 2 Se3, 
and Bi2Te3 [31 H], have been experimentally realized [5J- 
E]. The search for topological insulators has been ex- 
tended to ternary Heusler [9] and chalcogenide com- 
pounds [TO], but they are still limited to either the nar- 
row or zero gap semiconductors. Recently the layered 
honeycomb lattice Na 2 Ir03 [TJJ and pyrochlore A 2 Ir 2 07 
[El [13] have been suggested as possible topological in- 
sulators, although topological insulators with transition 
metal d electrons have not been fully investigated. Con- 
trary to the s-p electron systems, transition metal ox- 
ides with localized d electrons are expected to have both 
strong on-site Coulomb interactions and spin-orbit cou- 
plings. In particular, 5c? transition metal oxides such 
as iridates have a relatively weaker Coulomb correlation 
competing with spin-orbit coupled band structures. 

The interplay between the nontrivial band topology 
and the electron correlation effect can be an interest- 
ing development in the study of topological insulators. 
Ir-based transition metal oxides have shown some no- 
ble physics, such as the j' e ff = 1/2 insulating state 
in Sr 2 Ir04 |14H16j . the anomalous metal-insulator tran- 
sition in A 2 Ir 2 07 [17) . and the spin-liquid state in 
Na 4 Ir 3 08 [18] . where both spin-orbit coupling and corre- 
lation play important roles. 

Na 2 IrC>3 [HI [2D] has been proposed as a layered quan- 
tum spin Hall (QSH) insulator. Assuming the j e g = 1/2 
character of the band near the Fermi level [TT] , a single- 
band tight-binding model for the Ir-0 layer is mapped 
to the Kane-Mele model [1] for the QSH effect. The pro- 
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posed model Hamiltonian is, however, inconsistent with 
the first-principles band structure result which leads to a 
different prediction on the band topology of Na 2 Ir03 [21j . 
The inconsistency is partly due to the structural param- 
eters used for the first-principles calculations. It implies 
that its band topology may be sensitive to the structural 
variation. Consequently a small change of structure or in- 
teraction strength can drive a quantum phase transition, 
e.g., the change of its topological character. To under- 
stand the basic physics determining the band topology of 
Na 2 Ir03, we need to clarify the topological character of 
the spin-orbit coupled ground state and its dependence 
on the structure and interaction strength. 

In this Letter, we present a quantum phase transition 
between topological insulators (TIs) and normal insula- 
tors (NIs) in Na 2 IrC>3 based on an effective tight-binding 
Hamiltonian. We derived the effective Hamiltonian from 
the realistic tight-binding description of first-principles 
band structures. The electronic structure of the pla- 
nar edge-shared Ir0 6 octahedra contains large trigonal 
crystal field and direct hopping terms as well as a signifi- 
cant long-range hopping between extended Ir 5d orbitals. 
The phase boundary between topological and normal in- 
sulating phases is shown to depend on both the trigo- 
nal crystal field and the long-range hopping in Na 2 Ir0 3 . 
From the first-principles calculations of model structures, 
which simulate the change of tight-binding parameters, 
we confirmed that the interlayer distance can play a cru- 
cial role in the determination of the three-dimensional 
strong TI in Na 2 Ir0 3 . 

We introduce a tight-binding (TB) model for Na 2 IrC>3 
from the results of first-principles calculations [22 . The 
model is based on the Ir i 2g manifold in the two- 
dimensional honeycomb lattice and incorporates param- 
eters for an indirect hopping through the oxygen 2p- 
orbital (t p d), two kinds of direct hopping between neigh- 
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FIG. 1. Hopping parameters considered in our tight-binding 
model : (a) the indirect hopping(t p d) mediated by the oxygen 
2p orbital and two kinds of direct hoppings (td dl and t dd2 ), 
(b) the second-nearest-neighbor hopping (t ni ) and the third- 
nearest-neighbor hopping {t„ 2 ), and (c) the band structure of 
first-principles (solid lines) and tight-binding (dashed lines) 
calculations with spin-orbit coupling. 



Parameters Character Value (eV) 

t pd (pd7r)7(e d - e p ) 0.25 

t ddl Udda) + UddS) -0.5 

t dd2 l(ddn) + \{dd5) 0.15 

t n (sda) 2 /(e d - e„) -0.075 

A E R , g - E aig 0.6 

A Spin-orbit coupling 0.5 



TABLE I. Values of tight-binding parameters (in eV) obtained 
from the first-principles band structure. 



boring Ir atoms {tddx an d tdd 2 )i ano - another indirect hop- 
ping through the sodium 3s-orbital (t n ), as illustrated 
in Fig. 1(a) and (b). It is noted that the indirect hop- 
ping through the sodium 3s-orbital of Fig. [ljb) makes 
the second- and third-nearest-neighbor hopping terms 
non-negligible and plays a crucial role in the determina- 
tion of the topology of the band structure. In addition, 
we take into account both trigonal crystal field A and 
spin-orbit coupling (SOC) A terms. The TB parameters 
are summarized in Table |TJ The energy unit is in eV. 
The values of the parameters are fitted to the results of 
density-functional-theory (DFT) calculations within the 
local density approximation (LDA). For the DFT calcu- 
lations, we used the DFT code, OpenMX [23], based on 
the linear-combination-of pseudo-atomic-orbitals method 
|24j . The SOCs were included via the relativistic j- 
dependent pseudopotential scheme in the noncollinear 
DFT formalism [2"5rf2"7] . The unit cell we used in the 
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TABLE II. Parity invariants S(Ti) and Z2 topological invari- 
ants (v) with and without t n as determined from the product 
of parity eigenvalues at each time-reversal-invariant mome- 
tum T, Mi, Af 2 , and M 3 . 



DFT calculations is a simplified version of the original 
crystal, where the c-axis periodicity is reduced by chang- 
ing the relative stacking of Na layers with respect to the 
Ir network and also by neglecting the distortion of oxygen 
octahedra. We have checked the band structures with the 
different stacking of Na layers and the rotation of oxy- 
gen atoms around the axis perpendicular to the plane, 
and observed only a small change in the band dispersion 
and fitting parameters. Our TB model describes well the 
DFT band structure, especially the bands near the Fermi 
level (Ep), as shown in Fig. [ijc). 

It is interesting that the calculated parity invariants 
and Z2 topological numbers for our TB model depend 
critically on the second- and third-nearest-neighbor hop- 
ping term t n , as demonstrated in Table |Tl) The topolog- 
ical invariants were determined by following the method 
proposed by Fu and Kane [28] . The Z2 topology of the 
DFT band structure with t n — —0.075 turns out to be 
trivial. This result is consistent with the previous first- 
principles calculations result [21] . contrary to the quan- 
tum spin Hall insulator phase predicted by Shitade et al. 
[11] , This discrepancy is likely due to an over simplifi- 
cation of the tight-binding model employed in Ref . [TT] . 
Their model was based on the assumption that the low- 
energy degrees of freedom are determined by the half- 
filled j e g=l/2 doublets. However, when the significant 
trigonal crystal field is introduced in Na2lr03, j c s = 1/2 
and j e ff= 3/2 are not well separated and the j c ff=l/2 
doublets no longer serve as a useful basis [21]. The ex- 
tended nature of the 5d-orbital combined with the edge- 
shared octahedral structure in the L^/sNa^C^ plane of 
Na2lrC>3 is a source of such a strong trigonal crystal field. 

The topological character of Na2lr03 is quite sensitive 
to the magnitude of t n . By turning off the second- and 
third- nearest-neighbor hopping, i.e., t n — 0, the system 
becomes a nontrivial TI with v = 1. The tiny differ- 
ence in the magnitude of t n is responsible for the inver- 
sion of the valance and conduction bands, thereby lead- 
ing to the change of the S(M) sign. It indicates that 
t n is a key control parameter for the 'band-inversion' in 
this system. Therefore, the TB model with the nearest- 
neighbor hoppings only can not be sufficient for the de- 
scription of the topological character of Na2lrC>3 even 
though the nearest-neighbor TB models can describe rea- 
sonably the band dispersions of the j c s = 1/2 states in 
Sr 2 Ir0 4 OUS] and the hyper-Kagome Na 4 Ir 3 8 [HHH]. 

The nontrivial (trivial) Z 2 topological number can also 
be confirmed by the odd (even) number of pairs of gapless 
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FIG. 2. One-dimensional energy bands for armchair strip in 
the (a) NI phase with t n = -0.075 and (b) TI phase with 
in = 0. 
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FIG. 3. TB band structure with the variation of t n . The 
solid (red) (t n = 0) lines belong to the TI, and the dot- 
dashed (blue) (t n = —0.040) and double-dashed (orange) 
(t n — —0.075) lines belong to the NI. Note that, at the tran- 
sition point, i.e., t n = —0.021, the dispersion is linear. 



edge states. To examine the edge states, we constructed 
a TB Hamiltonian for the strip geometry with two edges 
in an armchair configuration. Figure [2] shows the one- 
dimensional energy bands with t n = —0.075 and t n = 
0. The bulk states and gaps are clearly seen and there 
are edge states which transverse the gap. In the normal 
insulating phase (u = 0), the edge states cross the Fermi 
energy an even number of times, as expected from the 
trivial Zi topological number. 

Transitions between trivial and nontrivial phases can 
be tuned by adjusting the key parameters of the band 
structure. Murakami and co-workers |31| developed the 
low-energy effective theory for the phase transition be- 
tween TI and NI systems in 2D and discussed the clas- 
sification of the possible types of transition. Following 
their classification, the gap closing at the time-reversal- 
invariant-momenta (k = G/2) occurs in systems with 
inversion symmetry. In our model, t n controls the gap 
closing at the M point in the Brillouin zone contrary 
to the case of a honeycomb lattice model such as the 



Kane-Mele model, where the gap closes at the K and K' 
points. The low-energy long-wavelength effective Hamil- 
tonian, which is expanded to linear order in k around the 
M point, can be written as 
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where k' = k — M and z± = b\k! x + bzk' y ± ib 2 k' y with 
real constant b\, 62 and 63 [3T]. Since the eigenenergies 
E = E ± yjm 2 + z + z_, m — corresponds to the gap 
closing. 

Figure [3] shows the TB band structure along some high- 
symmetry line for the several values of t n . In the TI re- 
gion (t n = 0) , the valance and conduction bands are sep- 
arated by a finite gap with negative mass. When \t n \ in- 
creases, the direct gap at the M point decreases. Finally, 
the gap collapses at the transition point t n ~ —0.021. As 
\t n \ passes through the transition point, the gap reopens 
and the system becomes a normal insulator with posi- 
tive mass. During this procedure, the order of the bands 
is inverted around the M point, which characterizes the 
topological nature of the system. 

To further verify the relation between the mass and 
the control parameters, we obtain the expression for the 
mass as a function of t n , A and SOC (A): 



2m = A + BA 



DAt ni . 



(2) 



Here, to clarify the role of second- and third-nearest 
hoppings, we refined t n into the second-nearest-neighbor 
hopping, t ni , and the third-nearest-neighbor hopping, 
t n2 . The expansion coefficients A, B, C, and D can be 
represented as functions of A. In the range of A from 
O.8A0 to I.6A0, where Ao is the SOC determined by the 
LDA calculation, the coefficients can be fitted as follows: 



A = +0.0746 — 0.1438A 
B = -0.0237 - 0.3061A 
C = -2.2178 + 1.2660A 
D = +1.0975 - 0.6837A. 



0.2587A^ 



(3) 



In the limiting cases of X/Xq = 1 and 00, the coeffi- 
cients are also shown in Table III From Eq. |2]), we 
can determine the Z2 topologica 
for given t n , A, 



number of the system 
and A by checking the sign of mass. 
This mass function correctly characterizes the topologi- 
cal phase transition. 

Figure [4] shows a, X — t n phase diagram for several val- 
ues of A . For a sufficiently large A, an indirect bulk gap 
opens and the system becomes an insulator. The insu- 
lating region is divided into two phases: TI and NI. The 
phase boundary between TI and NI corresponds to the 
m = line. There seems to be an asymptotic limit of 
t n beyond which the system stays as a normal insulator 
regardless of the strength of the SOC A. The vertical 



4 



A/A 


A 


B 


c 


D 




S(T) 


8(M) 


6(A) 


5(L) 




1 


-0.00 


-0.12 


-1.61 


0.58 


c/c =1.0 


-1 


-1 


-1 


-1 


0; (000) 


oo 


-0.13 
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-1 
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-1 
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1; (000) 



TABLE III. The expansion coefficients of mass function in 
Eq. [2] when A/Ao = and 00. These coefficients at A/Ao = 1 
and IS. = 0.6 correctly describe the transition point t n ~ 
-0.021. 
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FIG. 4. Phase diagram as a function of t n and SOC for several 
values of A. The pink circle indicates where the reality exists. 
\t c \ = 0.062 is the asymptotic line in the case of A = 0.6. 



dashed line in Fig. [4] is described by 

.4 



tr. 
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-0.067 
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with A=0.6. This implies that an arbitrarily large SOC 
cannot guarantee the nontrivial topology even if the up- 
per four bands originate from the j e ff = 1/2 states. 

Another important factor controlling the topological 
character is the trigonal crystal field, A. In fact, both 
SOC and the trigonal crystal field are two competing 
parameters characterizing the bands near the Fermi level. 
While SOC prefers the formation of the j c g = 1/2 state, 
the large trigonal crystal field becomes an obstacle for 
the j e g = 1/2 state. In this system, however, the trigonal 
crystal field seems to play a crucial role. As far as B < 
and D > of Eq. ^ as shown in Table III it is clear that 
a large trigonal crystal field is in favor of the nontrivial 
topology. This trend in the parameter dependence of 
the topological character is clearly reflected in the phase 
diagram of Fig. [3] 

Returning to a realistic system of Na2lr03 materials, 
we tried to probe possible TI phases by performing first- 
principles calculations. According to the simulated TB 
parameters for the TI phase, we choose two represen- 
tative structures among the various structures: (i) the 
original geometry of Na2lr03 (c/c = 1.0) and (ii) a vir- 
tual structure with the interlayer distance expanded by 
30% (c/co = 1.3). Calculated Z 2 topological numbers 
for each configuration are as listed in Table [TV] Here the 
c/cq parameter represents the change of interlayer dis- 



TABLE IV. Parity invariants S(Fi) and Z2 topological in- 
variants calculated from first-principles calculations. We con- 
sidered two different structures: (i) the original geometry of 
Na2lrC>3 (c/cq = 1.0) and (ii) a virtual structure with the 
interlayer distance enlarged by 30% (c/co = 1.3). 



tance relative to the original one. It is remarkable to find 
a nontrivial Z2 number for the structure with c/cq = 1.3. 
In other words, the increased interlayer distance drives 
the NI into the TI. To understand the change of elec- 
tronic structure between c/cq = 1.0 and c/cq = 1.3, 
we constructed the maximally localized Wannier func- 
tion (MLWF) [51 [33J and obtained the MLWF effective 
Hamiltonian. It is found that the trigonal crystal field en- 
hanced by the increase of interlayer distance contributes 
to the change of Z2 character for Na2lr03. When we 
increase the interlayer distance, the energy level of the 
interplane Na 3s state lowers. The second- or higher- 
order hoppings through the unoccupied 3s state of inter- 
plane Na contribute to the on-site Hamiltonian matrix 
elements. Since these matrix elements give rise to the 
energy separation between e' g and a\ g states, it can be 
interpreted that the trigonal crystal field contributes to 
the band inversion at the M point. In fact the band 
inversion occurs only at three M points in the k z = 



plane. However, as shown in Table IV the parity invari- 
ants 5(Ti) for c/cq = 1.3 are quite different for the k z = 
and k z = tt planes. It indicates that the c-axis hopping 
is also important for the existence of the TI phase in the 
three-dimensional system of Na2lr03. 

In conclusion we provide a microscopic picture for the 
topological phase diagram in Na2lr03 and identify the 
key control parameters based on the effective Hamilto- 
nian analysis. We predict that the TI phase of Na2lr03 
can be realized by controlling the long-range hopping and 
trigonal crystal field terms. Our first-principles calcu- 
lations for the simulated Na2&03 model structure sug- 
gest that the interlayer distance can play a crucial role 
in the determination of a three-dimensional strong TI. 
In practice, we propose two ways of driving a transi- 
tion from NI to TI: (i) epitaxial strain by the substi- 
tution of Na by other elements such as Li or a film 
growing on a appropriate substrate or (ii) intercalation 
of some molecules for the increase of the interlayer dis- 
tance (e.g., is the water-intercalated Nao.35CoO2T.3H2O 
showing the superconductivity 34J). Indeed our first- 
principles calculations demonstrate that the topological 
insulator phase of Li 2 Ir03 can be achieved by 2% in- 
plane lattice strain [35] . In addition, other structure 
manipulation techniques developed for superlattice and 
heterostructure may be adopted to control the NI-to-TI 
transition and to design the topological-insulator-based 
devices. A recent experiment observed an antiferromag- 
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netic (AFM) insulating behavior in Na 2 Ir0 3 [MHO]- The 
AFM ordering in the TI phase of Na2ir03 can be a can- 
didate of a topological magnetic insulator [351 [37J or a 
topological Weyl semimctal [38 below ordering temper- 
ature (Tjv). Above Tjv, correlations could enhance the 
SOC effects due to the suppression of effective bandwidth 
to stabilize the TI phase [TJ] . Further there are proposals 
for the topological Mott insulator having gapless surface 
spin-only excitations [HJIH]. Na2lrC>3 is likely to provide 
a new playground to study the effect of the correlation 



in TI. 
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